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en . 

0^ ■ We introduce the running coupling constant of QCD in the high temper- 

I . ature phase, g^{T), through a renormahzation scheme where the dimen- 

T— I , sional reduction is optimal at the one-loop level. We then calculate the rel- 

^T \ evant scale parameter, Ay, which characterizes the running of (f'{T) with 

■^ I T, using the background field method in the static sector. It is found that 

P^' ^t/Ams = e(^-E+^/22V(4vr) ^ 0.148. We further verify that the coupling g^{T) 

X,' is also optimal for lattice perturbative calculations. Our result naturally ex- 

D . plains why the high temperature limit of QCD sets in at temperatures as low 

'T^ \ as a few times the critical temperature. In addition, our Ay agrees remark- 

. 5^ ' ably well with the scale parameter determined from the lattice measurement 

^ ■ of the spatial string tension of the SU(2) gauge theory at high T. 
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I. INTRODUCTION 

At high temperatures QCD is expected to undergo a partial dimensional reduction [|I],3, 
namely static correlations at distances larger than the thermal wavelength (1/T) can be 
reproduced by a three dimensional Lagrangian, where only the static modes of the original 
theory are present. This reduced Lagrangian can be computed perturbatively up to a specific 
order in the QCD running coupling constant. In fact, non-perturbative infrared phenomena 
(e.g. thermo-mass generation) prevent the complete reduction, i.e. reduction to all orders 
in the QCD running coupling, from taking place 0. Consequently, observables can be 
reproduced only up to corrections of a specific order, before non-perturbative physics begins 
to dominate. 

Even though a complete dimensional reduction is not possible, the partial dimensional 
reduction of QCD still provides a simplified physical picture. However, phenomenological 
applications of this picture depends crucially on how high is the temperature above which 
this picture begins to take place. Since we expect corrections to vanish with some power of 
the QCD coupling, and since at zero temperature the asymptotic freedom starts dominating 
QCD physics at typical scales of about 10 to 20 times Aj^, one might anticipate that the 
reduced theory should become effective only for T ^ ^ms- 

Contrary to this expectation, there are strong evidences [^]-P that the dimensional re- 
duction picture is already valid at temperatures as low as two or three times the critical 
temperature Tc (the deconfining transition in the pure Yang-Mills case or the chiral restora- 
tion in full QCD). Since T^ is numerically not very different from Aj^, and considering that 
the QCD coupling constant only runs logarithmically, it is very surprising to find that the 
high-T regime of QCD starts at temperatures this low. 

This apparent puzzle can be solved with the observation that the effective scale parameter 
for the reduced theory. At, is renormalized after integrating out the non-static modes and 
becomes drastically smaller than A^jg. In fact, the definition of a scale parameter that 
characterizes the approach to the dimensional reduction regime implies the definition of 
a suitable coupling constant, g'^{T), that yields a sensible perturbative expansion at high 
temperature, i.e. an expansion whose coefficients contain minimal contribution from non- 
static modes. 

In this paper, we use the background field approach to define and compute the relevant 
coupling constant, and hence the scale parameter At- More specifically, we calculate the 
one-loop effective action for the background field in the static sector, and define the renormal- 
ization scheme by requiring that dimensional reduction be optimal for this gauge-invariant 
quantity. Furthermore, we verify by an explicit computation that this same renormaliza- 
tion scheme is also optimal for lattice perturbative calculations at high T, and therefore it 
provides a natural scale also for lattice simulations. 

In section |T| we introduce the renormalization scheme that defines the scale parame- 
ter At within the background field approach. In section |T|, we apply this definition and 
calculate At/Aj^. First we perform the calculation for the SU(A^) gauge theory in the con- 
tinuum using dimensional regularization; the effect of light quarks is also considered. Then 
we repeat the calculation for the pure SU(iV) gauge theory on the lattice in the Wilson 
formulation. In section 0, we compare our result with numerical determinations of At 
from a lattice measurement of the spatial string tension at high T. Section ^ is reserved for 



the conclusions. Several technical points pertinent to the lattice perturbative calculation at 
high T are discussed in the Appendix. 

II. DIMENSIONAL REDUCTION AND OPTIMAL RENORMALIZATION 

SCHEME 

The standard SU(A^) Yang-Mills gauge theory reduces at the tree level to the three 
dimensional Yang-Mills theory with adjoint Higgs (0" = Qq) 

Ckb = -\f^F^ - \{D4r{D4r , (1) 

where F^^ = diQ'^ - djQI - g-it^'^Q\Q'j and {D4Y = d^'" ~ g^f^'^'Q^A''- The coupling 
gs is related to the four dimensional coupling through g^ = g^T. Since £rd is a super- 
renormalizable theory in three dimensions and there is no other dimensionful scale around, 
all the dynamical scales must be set by the coupling constant g^ = g^T. 

Of course, once loop corrections are included the reduced theory in Eq. (0) would acquire 
new vertices and the coupling constant g^ would depend on the original coupling g"^ in a 
more complicated way. For example, g^ would receive corrections, such as g^T and so 
on. However, due to the asymptotic freedom of QCD {g"^ ~ 1/lnT) we still expect that 
dynamical scales are set by g"^ ~ g'^T, provided the scale parameter is chosen in a suitable 
way. 

Therefore, we believe that the concept of dimensional reduction involves two equally 
important aspects. On one hand, there is the possibility of a simplified description by using 
a theory £rd with less degrees of freedom in lower dimensions. On the other hand, the 
evolution of the parameters of £rd as a function of temperature should be dictated by the 
original theory. The main concern of our present work is to determine this evolution, which 
in turn determines the temperature dependence of the relevant physical observables. 

A. Background Field Method in the Static Sector 

It is well known that the effective action calculated using the background field method |^ 
is gauge invariant for the background gauge field at T = 0. This gauge invariance guarantees 
that the coupling constant renormalization is related to the wavefunction renormalization 

— 1 /2 

of the background field through Zg = Zj^ . Hence, the calculation of the quadratic part 
of the effective action, i.e. the two-point function for the background field, is sufficient to 
renormalize the coupling |^. 

Moreover, to the leading order, there is no magnetic mass generation at finite T. There- 
fore, the one-loop effective action for the magnetic sector is invariant under time-independent 

— 1/2 

gauge transformations also at finite T, insuring that the relation Zg = Zj^ still holds for 
the static background field 

A«(r,a.) = 0, Al{T,x)=Nl{x). (2) 

The same conclusion can also be reached more formally by applying, for instance, the meth- 
ods of Ref. P to the background field of Eq. (Q). The residual gauge invariance in the static 



sector implies that, in order to compute the couphng constant renormahzation at finite T, 
we only need to compute the two-point function of the background field Af in the static 
sector. 

We can still use the zero temperature Feynman rules, as given for instance by Abbott [§ . 
The only difference in the calculation is that time-components of all momenta become dis- 
crete Matsubara frequencies {27mT), and the corresponding integrals become discrete sums. 

B. Subtraction Scale 

As exhaustively discussed by Landsman 0, the decoupling of the non-zero modes at high- 
T is maximal only in some specific renormahzation schemes, such as the BPHZ scheme. In 
the background field method we only need to fix one renormahzation condition: we demand 
that the two-point function for the background field in the low external momentum (relative 
to T) limit coincides with the contribution solely from zero modes. 

Landsman [0] uses a finite temperature renormahzation group approach, since he dis- 
cusses thermal reduction in a more general context where several couplings are present. 
Thanks also to the background field approach, we deal with a simpler situation where only 
one coupling needs renormahzation. 

Therefore, we can directly implement the renormahzation condition by using the freedom 
in the choice of the subtraction scale, /i, which becomes a function of T. Intuitively, we expect 
fi to be of order of T. The purpose of our paper is to find out what is the proportionality 
constant. 

Then the reduced theory, Eq. (|T]), with the T-dependent coupling f^f = g^{fi{T))T, 
reproduces the full two-point function up to corrections of order of p'^/T'^ at the one-loop 
level. Due to the gauge invariance, the two-point function for the static background field 
A°^ must have the form 

(Sijp'^ - piPj ) Sab IIm (p^ , r, /i) • (3) 

Specifically, we choose fi by requiring the following renormahzation condition for the non- 
static contribution to IIm{p^,T, ^): 

UZ'{p' = 0,T,^^iT)) = 0. (4) 

The procedure is best explained by directly going through the calculation in the next section. 

III. CALCULATION OF THE SCALE PARAMETER 

A. In the Continuum 

In the continuum calculation we use dimensional regularization in the spatial dimensions, 
that is 



and the MS subtraction scheme. 

At the one-loop level there are four graphs that contribute to Um in the full theory: the 
bubble and tadpole graphs for both the quantum gauge fields and the ghost fields. The 
resulting total contribution is 
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where g'^{fi) is the running coupling defined in the MS scheme, f3o = llA^/(48'7r^), a is the 
gauge parameter, and 7^ is the Euler constant. 

The first term in Eq. (^) is obviously the classical contribution. The second term is the 
contribution of the static modes, and one can easily check that it can be reproduced by the 
reduced theory, Eq. (|I|), with coupling constant g^ = g^{^)T and the same gauge param- 
eter. The third term is the one that must be eliminated according to our renormalization 
prescription, which accomplishes maximal decoupling |0. We obtain this result by choosing 



/x(T) = AnTe 
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where c = 1/22. It is very reassuring to find that this optimal choice of the subtraction 
scale /i is independent of the gauge parameter a. The remaining contributions in Eq. (^ 
are suppressed by powers of 1/T^. 

In summary, to achieve maximal decoupling and hence the optimal dimensional reduc- 
tion, the effective coupling in the reduced theory must be 
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which defines the scale parameter 
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This results has a clear physical interpretation. The non-static modes decouple in the 
high-T limit, but their presence is nevertheless revealed by the appearance of the new scale 
At in the reduced theory (without any reference to the original theory, the only scale would 
be T). While this new scale is obviously related to the scale Aj^ that governs the full 
theory at zero temperature, the two scales do not coincide. Only with the coupling constant 
of Eq. (H), whose temperature evolution is set by the scale At in Eq. (||), the reduced 
theory is capable of reproducing the full four- dimensional one-loop corrections up to terms 
suppressed by 1/T^. Incidentally, it is interesting to notice that, if one intuitively identifies 
27iT (rather than T) as the relevant frequency unit, one gets an answer numerically close to 
the right one in Eq. (|^). 

At this point, we must point out that the optimal perturbative dimensional reduction 
criterion alone does not uniquely determine the scale A^. In general. At also depends on 
the specific Green's functions for which we demand optimal reduction. The use of a different 



process, represented by a different set of Feynman graphs, yields a different c in Eq. (^. 
However, we believe that typically \c\ ^ 1, and a different choice should not modify the scale 
ratio in Eq. (^) in an essential way. For example. Landsman calculated the temperature 
dependent coupling renormalization factor Zg by imposing maximal dimensional reduction 
on the two- and three-point functions in the conventional effective action (where the rela- 
tionship Zg = Z^ no longer holds). He did not express his result explicitly in terms of 
the scale ratio. But if we do it, we find that his result is quite close to ours, i.e. Eq. (|]) with 
c = 0. 

Another example that clearly shows the necessity of using an optimal dimensional re- 
duction scheme for defining the relevant scale at high-T can be found in the Gross-Neveu 
model . In that model a similar strategy makes the sub-leading correction to the screening 
mass of order of g^{T), rather than g'^{T), demonstrating that g^{T) is a sensible expansion 
parameter. 

Of course, the optimal dimensional reduction criterion is not the only way to define a 
temperature dependent coupling constant. For example, the quark- ant iquark potential at a 
distance of order of 1/T is used to define g'^{T) in Ref. 0. While it is certainly legitimate 
to make such a choice, it is also true that, because the reduced theory is meant to reproduce 
the full theory only at distances much larger than 1/T (spatial momenta small compared to 
T), definitions of the couplings made by matching short distance properties of the full and 
reduced theories do not necessarily define a scale that correctly characterizes the approach 
to the asymptotic high-T regime. 

At last, let us consider the effect of quarks on our result. If Nf light quarks are present 
in the theory, results of Eqs. (0), (D and (D still apply, but with /3o = (HA^ - 2Nf)/{A87r^) 
and c = {N/2 — 2Nf In4)/(11A^ — 2Nf), where we have adopted the convention for the trace 
of the Dirac- matrix: Trj^ji, = —46^^. 

For the phenomenological relevant case of A^ = Nf = 3, we get the value c = —0.2525, 
which corresponds to Ay/Aj^ ^ 0.110. Therefore, additional flavors further decrease the 
scale ratio until Nf > 16 (for A^ = 3), where asymptotic freedom is lost. 

B. On the Lattice 

The determination of the scale parameter A^ that governs the temperature dependence 
of the coupling in the reduced theory only involves an one-loop calculation. Nevertheless, 
the reduced theory is in general still non-perturbative, and non-perturbative methods are 
necessary to extract information from it. The standard approach is of course the lattice 
formulation. 

It should be clear from its definition (see section |T| and Ref. 0) that the concept of 
maximal decoupling scheme, along with the associated scale parameter, is independent of 
how the theory is regularized, and we expect the same Ay on the lattice, as long as we use 
the same renormalization condition, Eq. (^. 

On the other hand, since the lattice theory is usually defined in terms of the bare lattice 
coupling, Qoia), without introducing additional renormalization scale other than the lattice 
constant a, it is quite interesting to see with an explicit calculation how this scale emerges on 
the lattice at high T. In this respect, there are close analogies between our choice of g'^{T) as 
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a suitable expansion parameter for lattice perturbative calculation at high temperature, and 
the necessity of using expansion parameters different from the bare lattice coupling g^ia) 
for perturbative calculations at zero temperature p. 

In the following we verify that optimal dimensional reduction for the lattice effective 
action computed in the background field method defines indeed the same scale parameter 
we have found in the continuum calculation. For the sake of concreteness, we perform the 
calculation for the pure SU(A^) Wilson action, but the same result is expected to hold for 
other actions as well. 

In general, the coupling defined in the lattice background field method should have the 
following dependence on the bare lattice coupling up to one-loop 
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We want to show that c^ is such that g1{T) = g'^{T). Since we have expressed g'^{T) in 
terms of g'^{fi) in the MS scheme, see Eq. (H), we use the known relation between gQ{a) and 
g^in) in the MS scheme ||,|10| 
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and express also g^iT) in Eq. ( p!OD in terms of g^{fi) 
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By comparing Eq. (H) and Eq. (0) we see that to show g\{T) = g'^{T) is equivalent to show 
that 
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where the explicit expression of c° , which has been calculated by several authors P,pi]|, is 



c° = — ( -ll^E + 2/n + 3/oo + 6/io - 1 + 2Ati^z^^ + Gtt^ - Qti^/N^ 



The constants fij and Zij are defined as 
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Here and in the following we have closely followed the notation of Ref. 
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Since most of the calculation of the finite constant cj^ is closely parallel to the calculation 
of c^^, we only report the final result. In the Appendix, however, we illustrate the only new 
ingredient that is not a trivial extension of the calculation of c°: the high temperature 



expansion on the lattice. The lattice correspondent of the continuum result of Eq. 
Feynman gauge (a = 1) is 



^t,(p^T,a) 
with c^ given by 
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Since, as shown in the Appendix, flj = fij — 7^; — 3 In 4, this complete the proof of Eq. ([13D 
and, therefore, of the fact that g\{T) = g^(T). 

In other words, if we use g'^(T) as the expansion parameter, the lattice effective action 
in the high-T limit takes the following form 
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which is the same as its continuum counterpart, if we use the same coupling g'^{T) (see 
Eq. (^) with a = 1 and /i given by Eq. (0) ). In both cases we have been able to absorb 
in the coupling constant all leading local corrections due to non-static modes, while the 
non-local ones are reproduced by the reduced theory. 

IV. COMPARISON TO LATTICE RESULT 



In the preceding section, we have demonstrated that A^ is the relevant scale parameter 
in the high-T limit. Our argument is yet only perturbative in nature. However, as we 
emphasized earlier, the determination of the scale parameter is largely one-loop effect. Now 
let us compare our result with the scale parameter determined from a non-perturbative 
method: lattice measurement of the spatial string tension at high T. The primary reason 
for choosing the spatial string tension [^ rather than the heavy quark potential at distances 
of order of 1/T is that the concept of dimensional reduction only makes sense for large 
distance (low momentum) quantities. 

Bali et al. measured the spatial string tension in SU(2) gauge theory as a function of 
temperature o"s(T). Then they fitted their result to the expected form of the string tension 
in the three-dimensional SU(2) Yang-Mills theory 



a,{T)^g\T)T, 
where the running of g'^{T) with temperature is determined by the SU(2) /9-function 
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Even though the simulation process knows nothing about the dimensional reduction, the 
fitting formula Eq. (^) in fact defines the optimal three-dimensional coupling g^ = g'^{T)T 
through the string tension, similar in spirit to what we have done for the background field 
effective action. As a result, their fitted value of A^ = (0.076 ± 0.013)rc is the first, to our 
knowledge, non-perturbative determination of the scale that characterizes the high-T regime 
for the SU(2) gauge theory. 

In the scaling regime we expect that the critical temperature behaves like 



51 
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where jSc = 2N/gQ{a). From their critical coupling jSc = 2.74 at N^ = 16, and the known 
ratio Aj^/Al = 38.85 exp[—37r^/(llA^^)], it is straightforward to express Tc in terms of Aj^: 
Tc = 1.62A^jg. Then their numerical measurement yields 

AJ = (0.123 ±0.021)Aj^, (24) 

which is remarkably close to our result 

7E+1/22 

A^ = ^^— Aj^^0.148Aj^, (25) 

in spite of the different renormalization conditions. 

The smallness of A^/Tc gives a natural explanation of why their spatial string tension 
already at temperatures around 2Tc is numerically so close to the string tension of the three 
dimensional SU(2) Yang- Mills gauge theory. 

The result of Ref. @] implies a minor role for the Higgs sector in the reduced theory, 
whereas the result of Ref. [Q] seemingly implies the contrary. It would be interesting to see 
whether the use of an optimal coupling in calculations such as those in Ref. [^ and focusing 
only on long distance quantities could resolve the disagreement. For example, the gluonic 
Debye-screening mass, fin ~ g(T)T, could be used to define yet another optimal coupling. 
It is rather unfortunate that the numerical results in Refs. |]^JTT| are not accurate enough to 
determine a meaningful At- 

V. CONCLUSIONS 

We have defined a temperature dependent running coupling constant g^{T) in the spirit of 
the maximal decoupling of non-static modes of Landsman for SU(A^) gauge theories. More 
specifically, this coupling is such that the static modes reproduce the quadratic part of the 
one-loop effective action for the background field in the low momentum limit. Furthermore, 
g'^{T) provides a meaningful expansion parameter in the high-T limit, and its dependence 
on temperature defines a typical scale At for the high-T regime. 

We have calculated this coupling constant, and the related scale parameter, first in the 
continuum with dimensional regularization, where we verified its independence from the 
gauge fixing parameter a, and then showed that the same coupling is also optimal for the 



lattice perturbative calculation at high T. Our results are At = 0.148A|^ in the pure 
Yang-Mills or quenched cases and At = O.llOAj^ ioi N = Nf = 3. 

We have argued that this scale is typical in the high-T regime, even if its precise value 
depends on the specific definition. The consequence of our result is that the high-T regime of 
QCD, where the dimensional reduction picture appears to take place, sets in at temperatures 
as low as a few times of the critical temperature. 

Our calculation is in very good agreement with the non-perturbative determination of 
the scale parameter in the lattice simulations ||^ in the SU(2) Yang-Mills theory, therefore 
reinforcing the advantage of the renormalization scheme based on the optimal dimensional 
reduction criterion. 

It would be of great interest to have other lattice measurements of the scale parameter 
using other observables, such as the ones related to the gluonic Debye-screening mass and 
the deviations of the mesonic and baryonic screening masses from their free values. 

This work was supported in part by funds provided by the U.S. Department of Energy 
(DOE) under contract number DE-FG06-88ER40427 and cooperative agreement DE-FC02- 
94ER40818. 

APPENDIX: 

In this appendix we discuss several points of the high temperature expansion in per- 
turbative lattice calculations. First we use the ghost bubble-graph to illustrate the general 
method, then we prove that /^ — fij = —•je — 3 In 4, and finally discuss the convergence of 
the frequency sums to the corresponding zero temperature integrals. 

From the lattice action, see for instance Ref . 0] , we derive the following expression for 
the ghost bubble-graph 

JY (Q-ikfj,a _ ^i{k^-pfj.)aW^-i{k^-p^)a _ ^ik„a\ 

^'"^^^ ^ 4^ ? Ea(1 - cos k^a) Ep(l - cos(A;, - p,)a) ' ^^^~^ 

where Vt is the space-time volume and p = (0,p). Exponentiating the denominator and 
converting the spatial momentum sums into integrals (we work in the infinite spatial volume 
limit), we obtain 

AT A^T — 1 „7r rl-ih coo o 3 , 



--'^'-i;S^S/:|;^r«.e— .-^"^.nev- 



'T n=0 "-^ (2^) -^ A=l 

X I p~*('^M+'=i'-Pi''^) _|_ pi(kti+kv-Py.a) _ i{kfj,-ki,-p^a+pva) _ „~i{kii-k„)\ ('A2'l 

where 0a is imphcitly defined by tan 0a = Psin{pxa)/[a + (3 cos{pxa)]. Now we perform 
the spatial momentum integrals, yielding the modified Bessel functions. For the sake of 
concreteness, let us consider the component /^ = 1 and z/ = 2 



]\T Nr-l „oo ? r ^ 



T n=l 

i{<t>-i+<t>2-p2a) _|_ „i{(l>i+4>2~Pia) _ „i(0i-02-piii+p2a) _ -i(4>i-4>2) 
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In Eq. (^^) B' is just B without the n = term in the frequency sum. This static term is 
in fact the one that is directly reproduced by the reduced theory, and should be excluded 
from the contribution due to non-static modes. 

In the limit of \p\a <^ 1 and |p| <^ T (we are interested in the small lattice spacing and 
high-T limit), Eq. ([A3|) further simplifies 
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The expression in Eq. ( |X^ ) diverges in the limit A^,- -^ 00. We explicitly isolate its divergent 
part with the following subtraction 
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Now the first term is finite in the limit A^ — »• cxo, i.e. the limit a — >■ with aN^ = 1/T fixed, 
and it is equal to — A'piP2/ii/(487r^). We then use in the second term the expansion 



vr 



Nr-l 
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AT ^-^ ci 



Nr ^1 sin(7rn/Ar,) 



27ij + In 
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and obtain 
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Next we want to relate the finite integrals that are found in the lattice high-T expansion 
f-j, defined in Eq. (|1^), to the corresponding integrals that are found in the zero temperature 
calculation fij, defined in Eq. (|15D . Directly from their definitions, we find 



/■CO 

fij - fij = 1™ / d.x e 



e{x - 2) 



e ""loix) 



X 



X 



lim 



-Ei{-2e)-2Q_y2{l + e] 



(A8) 
(A9) 



The expression is finite, and we have introduced a convergence factor e~^^ in the integral 
only to be able to integrate separately the two terms. At last we obtain the desired result 
fij ~ fij ~ ~lE — 3 In 4 by using the small epsilon expansions of the exponential-integral 
function Ei{—2e) = 7^; + ln(2e) -|- C(e) and of the Legendre function of the second kind 
2Q-i/2(l + e) = - ln(2e) + 3 ln4 + 0(e). 

The last issue we would like to address in this appendix is the high-T expansion of those 
terms independent of external momentum, such as the tadpole graphs. Physically we expect 
that these terms cannot contain In(aT), since they do not contribute to the renormalization, 
and in fact they should eventually cancel out due to the gauge invariance or the lack of 
magnetic mass generation at the one-loop level. Therefore we should be able to factorize 
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any power dependence on T trivially, and take the continuum limit of the frequency sums. 
Mathematically, this is guaranteed by the fact that the convergence of the limit 

1 ^^~^ / 27rn\ /"i 
lim — V f cos = / dx f (cos27ix) (AlO) 

N^^oo Nr ^Q \ Nr J Jo ' ^ ' ^ ' 

is exponential, at least when f{z) can be expanded as a power series in z, which includes the 
cases we are concerned with. Note that the terms with n = should be included in these 
tadpole-like graphs, since they are not reproducible by the reduced theory. 
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